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ABSTRACT 

We describe constraints on a "standard" 8 parameter open cold dark matter (CDM) model from 
the most recent CMB and SNla data. Our parameters are the densities of CDM, baryons, vacuum 
energy and curvature, the reionization optical depth, and the normalization and tilt for both scalar and 
tensor fluctuations. We find that although the possibility of reionization and gravity waves substantially 
weakens the constraints on CDM and baryon density, tilt, Bubble constant and curvature, allowing e.g. 
a closed Universe, open models with vanishing cosmological constant are still strongly disfavored. 
Subject headings: cosmic microwave background — supernovae: general 



1. INTRODUCTION 

The currently most popular cosmological model has of 
order N — 10 free parameters. Upcoming CMB exper- 
iments hold the potential to measure these parameters 
with unprecedented accuracy (Jungman et al. 1996; Bond 
et al. 1997; Zaldarriaga et al. 1997; Efstathiou & Bond 
1998), especially when combined with galaxy redshift sur- 
veys (Eisenstein et al. 1998) and supernovae la (SN la) 
observations (White 1998; Tegmark et al. 1998). Bow- 
ever, these papers have also demonstrated the importance 
of fitting for all N parameters jointly, revealing subtle de- 
generacies by exploring the full A-dimensional parameter 
space. For this reason, there has been a persistent drive 
towards larger N when analyzing data. The first analy- 
ses based on COBE DMR used N = 2 parameters, the 
CMB quadrupole normalization Q and the scalar tilt n s 
of the power spectrum (e.g., Smoot et al. 1992; Gorski et 
al. 1994; Bond 1995; Bunn & Sugiyama 1995; Tegmark 
& Bunn 1995). Since then, many dozens of papers have 
extended this to incorporate more data and parameters, 
recent work including Bunn & White (1997); de Bernardis 
et al. (1997); Ratra et al. (1998); Hancock et al. (1998); 
Lesgourges et al. (1998); Bartlett et al. (1998); Webster et 
al. (1998); Lineweaver & Barbosa (1998ab); White (1998); 
Bond & Jaffe (1998); Gawiser & Silk (1998), and Con- 
taldi et al. (1998). The most ambitious analysis to date is 
that of Lineweaver (1998 - hereafter L98), jointly varying 
N = 6 parameters: n s , Q, the Hubble constant h and the 
relative densities ft c dm, and of CDM, baryons and 
vacuum energy. 

A realistic minimal cosmological model should include 
all physically well-motivated parameters. Yet even the 
heroic L98 analysis lacks three: gravity-wave (tensor) fluc- 
tuations, parametrized by a relative quadrupole normal- 
ization r and a tilt n t , and the optical depth r from reion- 
ization. In an inflationary context, gravity waves are just 
as natural as deviations from n s = 1, and we know that 
r > since the Universe was reionized before z — 5. It is 
therefore timely to extend this drive towards larger N by 
analyzing this "minimal" 9-parameter model space. That 
is the purpose of the present Letter. 



2. METHOD 

In principle, such an analysis is straightforward: com- 
pute the theoretical CMB power spectrum Ci with the 
CMBfast software (Seljak & Zaldarriaga 1996) at a fine 
grid of points in the A-dimensional parameter space and 
make x 2_n ts to the available power spectrum measure- 
ments in Figure 1. In practice, this is quite tedious. With 
M grid points in each dimension, M N power spectra must 
be computed. Lineweaver's impressive N — 6 analysis in- 
volved running CMBfast millions of times, corresponding 
to years of workstation CPU time, and with M ~ 20 as 
in L98, the amount of work grows by more than an order 
of magnitude for each additional parameter. Fortunately, 
the underlying physics (see e.g. Hu et al. 1997 for a review) 
allows numerical simplifications as described below. 

2.1. Parameter space 

We choose our 9-dimensional parameter vector to be 
P = (uJcdm,uJb,T,h,n k ,n s ,n t ,Q,r), where the physical 
densities io i = h 2 fli, i = cdm, b. The advantage of this 
parametrization (see Bond et a l. 1997; Eisenstein et al. 
1998) will become clear in |2.3| . flk is the spatial curva- 
ture, SO f2 A = l-fifc-ri c dm-^b = l-Qk-(u c dm+Ub)/h 2 . 

We choose our grid to cover the following parame- 
ter ranges: 0.02 < u cdm < 0.8, 0.003 < u b < 0.13, 
< r < 0.8, 0.2 < h < 1.3, < n k < 0.9, 0.5 < n s < 1.6, 
0.24 < n t < 1. This extends the L98 ranges some- 
what, since L98 reported high likelihoods near certain grid 
boundaries. To avoid prohibitively large M, we use a 
roughly logarithmic grid spacing for w C[ f m , ujb and h, a 
linear grid spacing for fl k , a hybrid for n s , n t , r and no 
grid at all for the normalization factors Q and r. 

Although a fairly fine grid is desirable for the likelihood 
analysis presented in Section 3, we find that we can at- 
tain sufficient accuracy by running CMBfast on a coarser 
grid and then interpolating the multipoles CV onto the fine 
grid. To prevent the resulting model file from exceeding 9 
gigabytes in size, we also use an adaptive mesh approach, 
complementing the global grid with a finer subgrid in the 
most favored regions of parameter space. 
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2.2. Separating scalars and tensors 

If we were to run CMBfast in the standard way, com- 
puting scalar and tensor fluctuations simultaneously, we 
would have to explore an 8-dimensional model grid since 
only Q drops out as an overall normalization factor. In- 
stead, we compute the scalar fluctuations Cf calar and the 
tensor fluctuations C t l ensor separately, normalize them to 
both have a quadrupole of unity, and compute the com- 
bined power spectrum as 

C e = Q 2 [Q calar + rC\ ensor \ . (1) 

We therefore only need to compute two 6-dimensional 
grids with CMBfast, one over {u) c dm,i w f>j T , h, f2fc, n s ) and 
the other over (uj cdm , uj bl r, h, flk,n t ). 

In addition, we reduce the dimensionality of our parame- 
ter space to 8 by imposing the consistency relation (Liddlc 
& Lyth 1992) 

r = -1n t . (2) 

It holds in all monomial inflation models satisfying the 
slow-roll conditions (the 2 nd relation n t = n s — 1 holds 
only for a small subclass). We do this merely because 
it is well-motivated and reduces error bars - it does not 
accelerate our calculations. 

2.3. Separating small and large scales 

The multipolc moments Ce for £ <C 100 correspond to 
fluctuations on scales outside the horizon at recombina- 
tion. This makes them almost independent of the causal 
microphysics that create the familiar acoustic peaks, i.e., 
independent of u> c d m and u>b- We therefore compute the 
power spectrum for £ < 100 with the fine grid restricted 
to (t, h, fifc, n s ) or (r, h, f^, nt), using only an ultra-course 
three-point grid for uj c( i m and u>b to pick up weak residual 
effects aliased down from larger I. We then fill in the rest 
of the uJcdm- and w^-values by interpolation. 

For the remaining (high £) part of the power spectrum, 
more radical simplifications can be made. First of all, the 
effect of reionization is merely an overall suppression of Ce 
by a constant factor e _2r on these small scales. Second, 
the effect of changing both Qk and h (and implicitly J7a) 
is merely to shift the power spectrum sideways. This is 
because the acoustic oscillations at z > 1000 depend only 
on ujcdm and u>b, and the geometric projection of these 
fixed length scales onto angular scales 9 in the sky obeys 
6(Cl m , £1a, h) oc l/d J 4(il m , £!a, zi ss ). Here dA is the lumi- 
nosity distance 



d A = {\ + z lss )-±—!-, k = y/\l-n m -QA\, (3) 
ti n 



Jo [(l + z') 2 (l + r> m z')-z'(2 + z')^A] 1/2 ' 

where S(x) = sinhx, x and sinx for open (Q m + < 1)> 
flat (Sl m + f^A = 1) and closed (fl m + f^A > 1) universes, 
respectively. We compute zi ss , the effective rcdshift of the 
last scattering surface, using the fit in Appendix E of Hu 
& Sugiyama (1996). 

f2 m and Vt\ also modify the late integrated Sachs- Wolfe 
effect, but this is important only for £ < 30 (Eisenstein et 
al. 1998). The only other effect is a small correction due to 



gravitational lensing (Metcalf & Silk 1998; Stompor & Ef- 
stathiou 1998), which we ignore here because of the large 
error bars on current small-scale data. To map the model 
(Q m ,Q* A ,h*) into the model (Q m ,Q^,h), one thus shifts 
its high £ power spectrum to the right by an ^-factor of 

6(n* m ,n* A ,h*)/e(n m ,n A ,h). 

We therefore adopt the following procedure. We com- 
pute the £ > 100 part of the power spectrum for a 3- 
dimensional grid over {aj cdm , u b , n s ) or (uj cdmi Ub, n t ). We 
extend this grid to include h and Qk by shifting it sideways 
as described, then merge it with the low £ grid by adjust- 
ing its normalization to match at i = 100. In addition 
to reducing the dimensionality of the grids computed with 
CMBfast, this approach has the advantage that only flat 
models need to be run for the high grid, with the (much 
slower) computations involving curvature and reionization 
only being required up to £ — 100. 

Extensive tests show that these approximations typi- 
cally reproduce the power spectrum to about 5% accuracy 
for generic models, i.e., substantially better than the cur- 
rent measurement errors. As data quality improves, the 
errors introduced by the above-mentioned approximation 
scheme can of course be continuously reduced to zero by 
refining the (w c <j m , Wf,)-grid for low £ and shifting the splic- 
ing point upwards from £ — 100. 
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Fig. 1. — Three models (see text) are shown together with the 
37 CMB data points, with the line weight and shading emphasizing 
those with small error bars. 

2.4. Data and likelihoods 

We use the compilation of CMB data and window func- 
tions of L98 with the addition of the new QMAP results 
(Devlin et al. 1998; Herbig et al. 1998; de Oliveira-Costa 
et al. 1998), from which we use the two points combin- 
ing both flights. The 37 band powers are shown in Figure 
1. We compute our likelihood function as L(p) oc e _x / 2 , 
where the \ 2 fit of the data to Ce(p) is computed as in 
L98. This procedure has a number of deficiencies as we will 
now describe. The probability distributions for the mea- 
surements are not Gaussian. In addition, the error bars 
for all experiments include a sample variance term which 
depends on p, and this dependence is rarely included ex- 
plicitly in quoted measurement results. A better (offset 
log-normal) approximation for the band-power likelihood 
is given by Bond et al. (1998), but for most experiments, 
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the additional parameter that it requires has unfortunately 
not been computed yet. Finally, if the likelihood function 
L(p) is a multivariate Gaussian, then one can show that 
marginalizing (integrating) over a subset of parameters is 
equivalent (up to an irrelevant normalization factor) to 
maximizing over them. We will follow L98 in doing the 
latter, since it is both simpler and avoids the unpleas- 
ant ambiguities of choosing a Bayesian prior — alas, with 
a uniform prior, our 9-dimcnsional normalization integral 
would not even converge. As we will see, our L is in fact 
highly non-Gaussian in some directions, which means that 
our confidence limits must be taken with a grain of salt. 
They also depe nd on the choice of Bayesian prior, as de- 
scribed in §3.2. 



However, to put these statistical issues in perspective, 
this author feels that an even more pressing challenge will 
be to test the data sets for systematic errors, e.g., by com- 
paring them pairwise where they overlap in sky coverage 
and angular resolution (Knox et al. 1998; Tegmark 1998). 

3. RESULTS AND CONCLUSIONS 
3.1. Best fit 

The best fit model is shown in Figure 1, and gives 
X 2 — 22.9. The probability of obtaining such a low x 2 - 
value with 37 — 8 = 29 effective degrees of freedom is about 
22%, so although CMB experimentalists have occasionally 
been accused of underestimating their error bars, we are 
closer to the opposite situation here. 

It is noteworthy that despite our large parameter space, 
the best fit model r = r = 0, ttk = 0.3, h 2 £l c d m = 0.2, 
h 2 flb — 0.025, h — 0.5 and n s — 1.0 (solid line in Figure 
1) is comparatively boring, preferring neither reionization, 
gravity waves nor tilt and giving rather conventional values 
of h 2 fli, and h. Much more exotic models are also allowed, 
however. If we restrict the parameter search to flat models 
(fi fc = 0), the best fit is r = 0.5, r = 0, h 2 Vl cdm = 0.35, 
h 2 Q, b = 0.04, h = 0.4 and n s = 1.4, dashed in Figure 1, 
where the high acoustic peaks that would be caused by the 
strong blue-tilting and the high baryon density are tem- 
pered by very early reionization. If we restrict ourselves 
to inflationary "vanilla" models with r = 0,^ = r = and 
n s = 1, the best fit is h 2 fl cdm = 0.1, h 2 fl b = 0.02 and 
h = 0.4, dotted in Figure 1. 



3.2. Single-parameter constraints 

Constraints on individual parameters are shown in Fig- 
ure 2 and Table 1, interpolating their marginal distribu- 
tions. Gravity waves are seen to be generally disfavored, 
with the maximum-likelihood value n t = corresponding 
to r = 0, no gravity waves at all. The best fitting models 
all fail to quite match the low COBE DMR quadrupole, 
and tensors merely make this worse by adding additional 
large scale power. Reionization is also mildly disfavored, 
for the same reason — increasing r and simultaneously 
increasing Q by a factor e T causes mainly a net rise at 
small I. However, this feature is softer than that of grav- 
ity waves, so as illustrated in Figure 1, it can be largely 
offset by increasing n s , Lo b and ui cd m- The result is that 
there are no relevant constraints on r: not even the ex- 
treme case t = 0.8 can be ruled out from our CMB data. 

The thin lines show the constraints assuming r = r = 0, 
as in L98, and agree well with the L98 results consider- 



ing that these did not include QMAP. However, the heavy 
lines show that including r and r substantially weakens 
these bounds. Gravity waves and reionization soften the 
upper limits on n s , Lo c dm and u>b since they can lower the 
acoustic peaks given COBE-normalization on large scales. 
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Fig. 2. — Heavy lines show likelihoods for individual param- 
eters marginalized over all others. Thin lines show the stronger 
constraints resulting from assuming neither reionization nor gravity 
waves. If the likelihood were Gaussian, the 68% and 95% confi- 
dence limits would lie where the curves cross the two dashed lines 
(see Table 1). 

We have followed L98 in using a uniform prior, trun- 
cated outside the explored parameter range. Whereas L98 
limited this range to values considered reasonable, we have 
attempted to quantify what CMB alone can say, extending 
the range far enough for the likelihood to become small. 
Figure 2 shows that this was achieved for all parameters 
except uJcdm, where our exclusion of (quite unreasonable) 
values h 2 Q, c d m > 0.8 matters. For a full Bayesian analysis, 
our CMB likelihood function should be multiplied by the 
likelihood functions from all other relevant astrophysical 
measurements . 

3.3. Constraints on the acceleration of the Universe 

The above-mentioned fact that f2^ and h (and implic- 
itly J7a) both shift the high £ power spectrum sideways 
make them rather degenerate. To better understand the 
constraints on these quantities, we therefore plot them 
in the two-dimensional Q m — J7a plane (Figure 3), where 
Q m = £l c dm + Qb- Our results for t = r = agree well 
with those of L98 when considering that (a) our analysis 
includes QMAP and (b) we have plotted our 68% and 95% 
confidence contours at Ax 2 = 2.29 and 6.18, respectively, 
since they are two-dimensional, as in Press et al. (1992) 
§15.6, whereas L98 used A% 2 = 1 and 4. Unfortunately, 
CMBfast cannot currently handle closed (Qfc < 0) models 
(White & Scott 1996). As L98 points out, the likelihood is 
already decreasing as one approaches the diagonal fij, = 
line (dotted) from the lower left, so we have simply ex- 
tended our likelihood function to Qk < by extrapolation. 
When dropping the r = r = assumption, however, this 
is no longer true, and the upper right (light grey) region 
of the Q m — fl\ plane is no longer excluded. 

Figure 3 also shows that the constraints at the lower 
left are unaffected by reionization and gravity waves. This 
asymmetry is easy to understand physically. This region 
is ruled out because the first acoustic peak is too far to 
the right, whereas the light grey region had the peak too 
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far to the left. Adding a strong blue-tilt can shift the peak 
slightly to the right, but never to the left. Figure 1 showed 
that such a tilted peak could be lowered back to the origi- 
nal height using r and r, making it fit the data, but r and 
r clearly cannot raise a (red-tilted) peak. 
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Fig. 3. — Constraints in the £7 m — plane. Apart from our 
CMB constraints, Q m < 0.1 would be inconsistent with the amount 
of matter observed dynamically. For the age we have simply taken 
Hotg > 0.6 as a lower limit. The region preferred by SN la is 
that computed by White (1998) from the combined data of the two 
supernova teams. In addition, gravitational lensing constrains the 
upper left corner. 



The recent constraints from SN la are highly comple- 
mentary to our CMB constraints. Figure 3 shows the 
SN la constraints computed by White (1998) in a joint 
analysis of the published data from the two search teams 
(Perlmutter et al. 1998; Riess et al. 1998; Garnavich et al. 
1998). We see that even including r and r, the combined 
CMB and SN la data prefer Oa > 0.5, with a vanishing 



cosmological constant strongly disfavored. The conclusion 
that f2 TO < 0.5 does not survive the inclusion of r and r, 
however, and we cannot rule out the possibility that the 
Universe is closed. 

3.4. Outlook 

In conclusion, we have performed a brute force 8 param- 
eter fit of cosmological models to the currently available 
CMB data and compared this with SN la constraints. We 
found that although the inclusion of reionization and grav- 
ity waves weakened many bounds, interesting constraints 
remain on e.g. Quoted error bars on parameters have 
grown steadily since the first COBE results, as more pa- 
rameters have been included in the analysis. Since we have 
now extended our parameter space to essentially the full 
"minimal cosmological model" , the error bars might be as 
large in this Letter as they will ever get. From now on, 
the rapid improvement in data quality will hopefully de- 
crease them faster than they arc diluted by the addition 
of further parameters, ushering us into the era of precision 
cosmology. 

Table 1 - Maximum-likelihood values and 68% confidence limits 
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